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ABSTRACT 

We show that the predictions of the recent physical model by Granato et al. (2001) 
are fully consistent with the available estimates of the clustering properties of galaxies 
at (sub)-mm wavelengths. The model is then exploited to compute the expected power 
spectrum of fluctuations due to clustering at the frequencies of the High Frequency 
Instrument (HFI) of ESA's Planck satellite. The clustering signal is found to be 
detectable in regions of low interstellar dust emission. A further distinctive prediction 
of the adopted model is a remarkably high fraction of gravitationally lensed sources 
at bright mm/sub-mm fluxes. In fact, according to the adopted model, only thanks to 
gravitational lensing, a significant number of high-z forming spheroidal galaxies will be 
detectable by large area, shallow surveys at mm/sub-mm wavelengths, such as those 
carried out by Planck/HFI. Allowing for other source populations, we find that the 
fraction of gravitationally lensed sources at bright mm/sub-mm fluxes is expected to 
be up to ~ 40%. 
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1 INTRODUCTION 

The formation and early evolution of massive spheroidal 
galaxies is still a most controversial issue. In the general 
framework of the hierarchical clustering picture for structure 
formation in Cold Dark Matter (CDM) cosmologies, two sce- 
narios are confronting each other. According to the "mono- 
lithic" scenario, a large fraction of giant elliptical galaxies 
formed most of their stars in a single gigantic star-burst at 
high redshift, and then essentially underwent passive evolu- 
tion. This approach is in contrast with the "merging" sce- 
nario, backed up by semi- analytic galaxy formation mod- 
els (e.g. Kauffmann 1996; Baugh et al. 1998), wherein large 
galaxies mostly formed via mergers of smaller galaxies that, 
at the time of merging, had already converted part of their 
gas into stars. These two scenarios yield very different pre- 
dictions of the evolution of large ellipticals: if "monolithic" 
collapse is assumed, the comoving number density of such 
galaxies remains essentially constant with redshift, while 
their bolometric luminosity increases with look-back time 
out to their epoch of formation. On the contrary, accord- 



ing to the "merging" scenario, the comoving density of large 
ellipticals decreases with increasing redshift. 

The evolution of the comoving density in the optical 
band has however proven difficult to measure, and discor- 
dant results have been reported. A deficit of old ellipticals 
at z > 1, when compared to predictions from pure lumi- 
nosity evolution models, has been claimed by several groups 
(Kauffmann et al. 1996; Zepf 1997; Franceschini et al. 1998; 
Barger et al. 1999b; Menanteau et al. 1999; Treu & Stiavelli 
1999; Rodighiero et al. 2001), while others found evidence 
for a constant comoving density up to « ~ 1.5 (Totani & 
Yoshu 1997; Benitez al. 1999; Broadhurst & Bowens 2000; 
Scodeggio & Silva 2000; Cohen 2001). 

Furthermore, Daddi et al. (2000) find that the observed 
surface density of Extremely Red Objects (EROs)- the bulk 
of which appears to be passively evolving ellipticals - is con- 
sistent with predictions from the "monolithic" scenario. 

The results of SCUBA and MAMBO surveys at sub- 
mm/mm wavelengths severely challenge the "merging" sce- 
nario. The most exhaustive investigation carried out so far 
(Granato et al. 2000) combines the semi-analytic galaxy for- 
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mation model by the Durham group (Cole et aL 2000) with 
the state-of-the-art spectro-photomctric code by Silva et al. 
(1998) which includes reprocessing of starlight by dust. It 
turns out that, while this approach successfully accounts 
for a largo variety of obscrvablcs (spectral energy distribu- 
tions of galaxies of different morphological classes, luminos- 
ity functions over a wide wavelength range, from UV to far- 
IR, and more), it nevertheless falls short by a substantial 
factor when trying to account for the mm/sub-mm counts. 
Similarly, Dcvricndt & Guidcrdoni (2000), who used a to- 
tally independent semi-analytic model, were forced to intro- 
duce an ad-hoc galaxy population to reproduce such counts. 

The main problem stems from the fact that, although 
only a handful of redshifts of SCUBA/MAMBO sources have 
been reliably measured so far, there is growing circumstan- 
tial evidence that many (perhaps most) of them are ultra- 
luminous star-forming galaxies at a > 2 (Small et al. 2000; 
Dunlop 2000; CariUi et al. 2000; Bortoldi ct al. 2000; Scott 
et al. 2001; Fox et al. 2001). Furthermore, the inferred star- 
formation rates are very high (typically from a few hundreds 
to ~ 1000 M© yr~^), consistent with those required to build 
the stellar populations of massive ellipticals on a timescale 
< 1 Gyr. 

A detailed, astrophysically grounded, model which fully 
accounts for both the SCUBA/MAMBO counts in the 
framework of hierarchical clustering scenarios, and for the 
main aspects of the chemical evolution of spheroidal galax- 
ies (stellar mctalicity, luminosity-mctalicity rc^lationship, a 
enhancement) has been recently presented by Granato et al. 
(2001). Key ingredients of this model are: i) a tight relation- 
ship between formation and early evolution of spheroidal 
galaxies and active nuclei hosted at their centres (note that 
first evidences for a SCUBA/AGN connection have recently 
been found by Almaini et al. 2001); ii) once the effects of 
cooling and heating (the latter being mostly due to stel- 
lar feedback) arc properly taken into account, the timescale 
for star formation within virialised dark matter (DM) ha- 
los turns out to be relatively short for massive spheroidal 
galaxies while it is longer in the case of less massive halos, 
where the feedback from supernovae (and/or the active nu- 
cleus) slows down the star formation and can also expel a 
significant fraction of the gas. It therefore follows that the 
canonical scheme implied by the hierarchical CDM scenario 
- small clumps collapse first - is reversed in the case of 
baryons. 

During the active star-formation phase, largo spheroidal 
galaxies show up as luminous sub-mm sources, accounting 
for the SCUBA and MAMBO counts. The subsequent evolu- 
tion is essentially passive, except possibly for minor mergers, 
much like that expected in the "monolithic" scenario. 

As mentioned by Granato et al. (2001), the model im- 
phes that SCUBA/MAMBO galaxies are highly biased trac- 
ers of the matter distribution. The model also establishes 
well defined links between these galaxies and other galaxy 
populations such as Lyman Break Galaxies (LBGs) and Ex- 
tremely Red Objects (EROs). Studies of their clustering 
properties therefore provide a tost for it. A first discussion 
of this aspect was offered by Magliocchetti et al. (2001a, 
henceforth MA2001). In this paper we extend the analy- 
sis and compare the clustering properties predicted by the 
model with recent observational measurements at 850 /xm 
and 170 /Ltm. 



Another distinctive feature of this model is the ex- 
tremely steep, essentially exponential decline of the bright 
tail of mm/sub-mm counts of forming spheroids. This is due 
to the fact that, according to the adopted model, the bulk 
of star formation in such objects is essentially complete at 
« ~ 2. The bright counts therefore essentially reflect the high 
luminosity tail of the luminosity function. This in turn relies 
on the fact that within the hierarchical clustering scenario 
massive halos are exponentially rare at high redshifts. Hints 
of such a steep decline can possibly be discerned in the re- 
cent MAMBO (Carilh et al. 2000; Bertoldi et al. 2000) and 
SCUBA bright surveys (Scott et al. 2001; Borys et al. 2001). 

As extensively quantified and discussed in Sections 4 
and 5, such a steep decrease in the number counts for fluxes 
lOniJy < Sgsoyjm < 100 mjy implies a large fraction 
of strongly lensed galaxies to appear at bright mm/sub- 
mm fluxes. This fraction is found to be larger than those 
worked out for other (phenomenological) models that also 
successfully account for SCUBA/MAMBO counts (Rowan- 
Robinson 2001; Takeuchi et al. 2001; Pearson et al. 2001; 
Blain et al. 1999b). 

The layout of this paper is as follows. In § 2 we compare 
the predictions of the model by Granato ct al. (2001) with re- 
cent measurements of clustering properties of SCUBA galax- 
ies (Peacock et al. 2000; Scott et al. 2001) and draw some 
conclusions on the cell-to-cell fluctuations of the projected 
density of sources, found to be strong enough to account for 
the discrepancies in the number counts observed in different 
regions of the sky by different groups. In § 3 we present the 
power spectra of temperature fluctuations due to clustering 
of forming spheroidal galaxies in Planck/HFI channels and 
derive predictions for the power spectrum of extra-galactic 
background fluctuations at 170 /im, to be compared with the 
signal detected by Lagache & Puget (2000). § 4 gives the es- 
sential details of our calculations on the effect of lensing on 
source counts, while in § 5 we report our predictions for 
counts of lensed galaxies at mm/sub-mm wavelengths, with 
reference to those covered by the Planck/HFI instrument. 
Conclusions are presented in § 6. 

Throughout the paper we will assume a Cold Dark 
Matter flat cosmology with = 0.7, Ot/i^ = 0.025, 
Ho = 100ftokms~^Mpc'^ with ho = 0.7 and a COBE- 
normalised as = 1. We note however that our conclusions 
are only weakly dependent on the underlying cosmological 
model. 



2 CLUSTERING OF HIGH-Z FORMING 
SPHEROIDAL GALAXIES: MODEL 
PREDICTIONS VS. OBSERVATIONS 

A crucial test for a galaxy formation model is its ability 
to correctly describe the large-scale properties of the pop- 
ulations under investigation. In MA2001 we have presented 
predictions for the clustering of SCUBA sources and their 
contribution to the total background fluctuations in the sub- 
mm band, based on the model by Granato et al. (2001). The 
model has been proven to provide a good fit to the clustering 
properties of LBGs (Giavalisco et al. 1998), corresponding 
in our framework to the early evolutionary phases of low-to- 
intermediate mass spheroidal galaxies. 

In this Section we will extend the analysis of MA2001 
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and compare the results to the best up-to-date clustering 
measurements at 850 //m. We start off with the theoreti- 
cal expression for the angular two-point correlation function 
w{e): 



w{e) 



J N^{z) fe,\ (M^in, z) (dz/dx) C(r, z) dz du 



dz\ 



(1) 



(Peebles 1980), where ^(r, z) is the spatial mass-mass corre- 
lation function (obtained as in Peacock & Dodds 1996; see 
also Moscardini et al., 1998 and Magliocchetti et al. 2000), 



X is the comoving radial coordinate. 



|2U/2 



(for 

a flat universe and in the small angle approximation), and 
N{z) is the number of objects within the shell {z,z -\- dz). 
The relevant properties of SCUBA galaxies are included in 
the redshift distribution of sources N{z), and in the bias 
factor feeff (-Mmin, z) . 

The effective bias factor feoff (Afmin, z) of all the dark 
matter haloes with masses greater than some threshold 
Mmin is then obtained by integrating the quantity b{M, z) 
(whose expression has been taken from Jing 1998) - repre- 
senting the bias of individual haloes of mass M - weighted 
by the mass function nscuBA(Af, z) of SCUBA sources: 



CD 

+ 



1.4 



1.2 - 



0.8 



0.6 



T — I — I — I — I — I — I — r 



J I I I I I I UJ I I I I I I \ L 



100 200 300 
0[Arcsecs] 



400 



Figure 1. Angular correlation function of 5 > 5 mjy SCUBA 
sources. Model predictions are shown by the dashed curves (where 
the higher one is for M^^^i^/ M^^^^ = 100 and the lower one for 
-'^^halo/A^sph = 10), while data points represent the Scott et al. 
(2001) measurements. 



feoff (Mmin, z) 



J^^^^^ dM h(M, z) nscuBA(M, z) 
dM nscvBA{M,z) 



(2) 



with nscuBA{M, z) = n{M, z) Tb /th, where n(M, z) dM is 
the mass function of haloes (Press & Schechter, 1974; Sheth 
& Tormen 1999), Tb is the duration of the star- formation 
burst and is the life-time of the haloes in which these 
objects reside. 

The bias factor fenjr in Eq. (^ and the angular correla- 
tion function of Eq. (hi) have then been evaluated for three 
mass ranges: i) masses in the range Mgph ^ 10^ — W^'^Mq, 
duration of the star formation burst Tb ~ 2 Gyr, and 
typical 850 /im fluxes S < 1 mJy; ii) masses in the range 
Mgph ~ 10^°-10"Mq and Tb ~ 1 Gyr; iii) masses in the 
range M^ph > W^^Mq, Tb ~ 0.5 Gyr, dominating the 
counts at 850 /xm fluxes 5* > 5 — 10 mJy (note that by Alsph 
we denote the mass in stars at the present time). 

Furthermore, we considered two extreme values for the 
ratio between mass in stars and mass of the host dark halo - 
A^fiaio/Afsph = 100 and Mhaio/Aispii = 10 - which encompass 
the range covered by recent estimates for massive objects 
(McKay et al. 2001; Marinoni & Hudson 2001). Granato et 
al. (2001) found that this ratio is 20 — 30 for the massive 
objects. 

According to our model, sources are expected to be 
strongly clustered, with a clustering signal which increases 
for brighter sources and higher values of Mhaio/Mgph. From 
the above discussion it follows that measurements of w{8) 
are in principle able to discriminate amongst different mod- 
els for SCUBA galaxies and in particular to determine both 
their star-formation rate, via the amount of baryonic mass 
actively partaking of the process of star formation, and the 
duration of the star-formation burst. 

Tentative evidence for positive clustering of SCUBA 
galaxies on scales of l'-2' has been found by Scott et al. 
(2001) for a sample of 38 5* > 5 mJy sources observed in 
two regions of the sky corresponding to the ELAIS N2 and 



Lockman-Hole East fields. Although such measurements are 
dominated by noise due to small-number statistics, it is nev- 
ertheless interesting to note that - as illustrated by Figure ^ 
- our predictions show full consistency with the data (kindly 
provided by S. Scott). 

Peacock et al. (2000) analysed the contribution of clus- 
tering of unresolved (S < 2 mJy) SCUBA sources to the 
850 /im background fluctuations detected in the Hubble 
Deep Field (HDF) North. They found some evidence for 
clustering of the background source population, consistent 
with an angular correlation function of the form w{9) = 
{6/6o)~^ , with e — 0.8 and 6o = 1". Figure |^ presents our 
model predictions for the angular correlation function w{6) 
derived in the case of SCUBA galaxies with S < 2 mJy; the 
redshift distribution spans the range 0.7 < z < 6. 

Our results are fully consistent with the constraints on 
structure in the sub- millimetre background set by Peacock et 
al. (2000) and not far from their preferred clustering model 
(dotted line in Fig. |), if Mhaio/Afsph = 100. The shaUower 
slope predicted by our model at small angles (0 ^ 4 x 10~^ 
degrees) ultimately stems from the fact that at high z we 
enter the regime of linear growth of density perturbations. 
This entails a flattening of the slope of the two-point spatial 
correlation function ^(r, z) for z > 2. As a consequence, the 
angular correlation function w{9) obtained as in Eq. (|l|) by 
projecting (,{r,z) over a wide range of redshifts, up to 2 ~ 
6, will mirror this tendency to flatten out. Since a redshift 
range 1 < z < 6 is also indicated by the analysis of Peacock 
et al. (2000), the adoption of a flatter slope for 'w{6) in their 
clustering analysis might have been more appropriate. 

Strong clustering substantially enhances cell-to-cell 
fluctuations of the surface density of detected sources. The 
angular correlation function w{9) is in fact related to the sec- 
ond moment (variance) of the galaxy distribution function 
via the expression (Peebles 1980; Roche & Eales 1999) 
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Figure 2. Angular clustering of unresolved (S < 2 mjy) sources 
at 850 /im. The solid curve is obtained for Mjjaio/^sph = lOi 
while the dashed one represents the case Mjjaio/^sph = 100. The 
dotted line corresponds to w(6) = (6/1")""'* (Peacock et al. 
2000). 



N + I — 



(3) 



where iV - mean count of sources in the solid angle u - rep- 
resents the Poisson noise arising from the discrete nature of 
the objects, and the normalised variance Y? can be written 

as 




Figure 3. Predictions for the rms fluctuations relative to the 
mean number of detected sources per cell of area lu = 
at 850 /im. Dashed and dotted lines are for sources brighter than 
10 mJy or 1 mJy, respectively. Higher curves of each type cor- 
respond to Mijaio/A^sph = 100, lower ones to A/jjaio/A/sph = 10. 
The solid lines show the Poisson contributions for the two flux 
limits. The lines are drawn only for angular scales such that the 
mean number of source in the solid angle u) = 0-^ is > 1. 



clustering amplitude predicted for SCUBA sources can eas- 
ily account for the differences found in the number counts 
observed in different sky regions by different groups (see e.g. 
Small et al., 1997; Hughes et al., 1998; Blain et al., 1999b; 
Barger et al., 1999a; Eales et al., 2000; Scott et al., 2001). 



= y w{e) dtoi duo-i. (4) 

In the case of square cells with uj — Q x Q, Eq. (^ can be 
expressed as the two-dimensional integral 

E^(e) = eM dx I w{6) dy, (5) 
Jo Jo 

with 9 — \J +y'^- 

The quantity in Eq. (jH]) has been evaluated for the three 
cases of low-, intermediate- and high-mass SCUBA galaxies, 
introduced in this Section and based on the Granato et al. 
(2001) model, and for the usual two different values of the 
Afhaio/Afsph ratio. Results for the rms fluctuations relative 
to the mean number of detected sources per ceU, ^ly^/N 
[see Eq. (^)] are presented in Fig. ^, where dashed and dot- 
ted curves correspond, respectively, to sources brighter than 
10 and 1 mJy. We have adopted a mean surface density per 
square degree of TV = 180 for a flux limit of 10 mJy, (see 
Scott et al. 2001), and of iV = 1.3 x lO" for a flux limit of 
1 mJy (Granato et al. 2001). Fig. ^ clearly shows that the 
contributions due to clustering are generally larger than the 
Poisson contributions, shown by the two long-short dashed 
lines, and in some cases by a large factor. For surveys cov- 
ering small areas of the sky the cell-to-cell fluctuations are 
large (up to ~ 100%). We therefore conclude that the high 



3 MULTI- WAVELENGTH ANALYSIS OF THE 
CONTRIBUTION OF CLUSTERING TO 
BACKGROUND FLUCTUATIONS 

The strong clustering of high-z forming spheroidal galaxies 
adds an important contribution to background fluctuations. 
Measurements of such a contribution can be informative on 
both the nature and the properties of sources below the de- 
tection limit. On the other hand, these fluctuations may 
have a signiflcant impact on the detectability of Cosmic Mi- 
crowave Background (CMB) fluctuations. 

MA2001 (but see also Scott & White 1999 and Haiman 
& Knox 2000) estimated the power spectrum of clustering 
fluctuations at 850 /im. In the present Section we will ap- 
ply the same analysis to different wavelengths, correspond- 
ing to the effective frequencies of the channels of the High 
Frequency Instrument (HFI) of the Planck ESA mission: 
2100 ^im ^ 143 GHz; 1380 ^im 217 GHz; 550 ^m 545 
GHz; 350 /im 857 GHz. We will also work out predictions 
at A = 170/im, the wavelength probed by the FIRBACK 
deep survey (Dole et al. 2001). 

The counts of spheroidal galaxies are obtained from 
the model by Granato et al. (2001). The counts of spi- 
ral and of star-burst galaxies have been estimated by ex- 
trapolating to each considered wavelength the 60 /tm lo- 
cal luminosity functions of the two populations, derived by 
Saunders et al. (1990), and evolving them in luminosity as 
L{z) = L(0)(1 + z)" , with a = 3.5 in the case of starbursts 
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Figure 4. Integral number counts of galaxies predicted by the model of Granato et al. (2001) at A = 2100/im (top left-hand panel), 
A = 1380/im (top right-hand panel), A = 850/im (middle left-hand panel), A = 550^m (middle right-hand panel), A = 350/^m (bottom 
left-hand panel) and A = 175/im (bottom left-hand panel). The solid lines show the total counts while the dot-dashed, dashed, and dotted 
lines are respectively obtained for the three populations of forming spheroidals, spirals and star-burst galaxies. Data in the ISSOjum panel 
are derived from MAMBO observations (Bertoldi et al., 2001); a flux density ratio 5i2oo/Si380 = 1-5 has been adopted to extrapolate 
to A = 1.38 mm the MAMBO flux densities at A = 1.2 mm (according to the model, such ratio varies between 1.4 and 1.6, depending 
on the details of the dust temperature distribution). Data at 175/im are from Dole et al. (2001). Data at 850//m are from Blain et al. 
(1999b) (open circles), Hughes et al. (1998)(open triangles), Barger et al. (1999) (open squares), Eales et al. (2000) (filled circles), Scott 
et al. (2001) (filled triangles), Borys et al. (2001) (crosses) 
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Figure 5. Predicted power spectrum of temperature fluctuations STi = + 1)Ci/2tt (in units of K) as a function of tlie multipole I 

at 2100 fim (top-left panel), 1380 ^m (top-right panel), 550 fim (bottom-left panel) and 350 fim (bottom-right panel), central frequencies 
of the Planck/HFI channels. The dotted curves are for the 5cr flux detection limit of each channel. Higher curves of the same kind are 
for Miiaio/^sph=100, lower ones for Afijaio/^sph=10- The solid lines represent the power spectrum of primary CMB anisotropics as 
predicted by a standard ACDM cosmology (A = 0.7, Qq = 0.3, ho = 0.7). The upper long-dashed curves represent the contribution from 
Galactic dust emission averaged all over the sky, while the lower ones correspond to high galactic latitudes (|f)| > 80°) only. 



and a = for spirals. The extrapolation in frequency has 
been carried out adopting the observed spectral energy dis- 
tributions of NGC6946 and NGC6090, respectively in the 
case of spirals and star-burst galaxies. 

The contribution to the counts of forming spheroidal 
galaxies, when compared to the case for spirals and 
star-burst galaxies, is found to decrease with decreasing 
wavelength because of the increasingly less favourable K- 
corrections associated with the former population. Spirals 
and starbursts consequently become more and more im- 
portant at shorter and shorter wavelengths, and eventually 
dominate the counts for A < 200/im. This is illustrated by 



Fig. ^, where we show the contribution to the integral num- 
ber counts (represented by the solid curves) of the three 
populations of spheroidal galaxies (dot-dashed lines), spi- 
rals (dashed lines) and starbursts (dotted lines) at different 
wavelengths. 

Figure ^ shows that - under our assumptions - star- 
bursts and spirals constitute a negligible fraction of the total 
number of sources at mm/sub-mm wavelengths, the domi- 
nant contribution coming from high-z forming spheroidal 
galaxies (the same which - according to Granato et al. 
(2001) - show up as SCUBA sources at A = 850/im), at 
all but very bright fluxes (with a flux threshold which varies 
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Figure 6. As in Figure ^ but for the 850 ^lui case. Dotted lines 
are for a detection limit Sd = 100 mjy while dashed ones for Sd = 
10 mJy. Higher curves of the same kind are for Miiaio/-'^sph=100, 
lower ones for Afiialo/-'\^sph=10. 



with wavelength) where, as aheady discussed in the intro- 
duction and extensively discussed in Section 4, the number 
counts of spheroidal galaxies experience an exponential de- 
cline. The situation is reversed at 175 /im (bottom right 
panel of Figure where at all but very faint fluxes, spi- 
rals and star-bursts dominate the predicted number counts. 
In the following analysis we will therefore only consider 
the population of spheroidal galaxies when dealing with the 
mm/sub-mm wavelength range. 

The angular correlation function of intensity fluctua- 
tions due to inhomogeneities in the space distribution of un- 
resolved sources (i.e. with fluxes fainter than some threshold 
Sd), C{6), can be expressed as the sum of two terms Cp and 
Cc, the first one due to Poisson noise (fluctuations given by 
randomly distributed objects), and the second one owing to 
source clustering (see MA2001 and De Zotti et al. 1996 for a 
detailed discussion). In the case of highly clustered sources 
such as forming spheroidal galaxies, the Poisson term Cp is 
found (see e.g. Scott & White 1999) to be negligible when 
compared to Cc- We can therefore safely assume C ~ Cc, 
whose expression is given by 



(6) 



with &cff (Mniin, z) defined by Eq. (0), x being the comov- 
ing radial coordinate, 5z = c/Hq u {u is introduced in Eq. 
(1) and c is the speed of light), Zmax the redshift when the 
sources begin to shine and z{Sd,L) the redshift at which 
a source of luminosity L is seen with a fiux equal to the 
detection limit Sd- The effective volume emissivity j'cff is 



expressed as: 

'"'<^{L,z) K{L,z) LdXogL, (7) 

"1>(L, z) being the luminosity function (per unit logarithmic 
luminosity interval), K{L, z) the K-correction and Lmax and 
L-mi-n respectively maximum and minimum local luminosity 
of the sources. 

Cc{d) in Eq. (^ has been evaluated separately for the 
three cases of low-, intermediate- and high-mass objects, 
by plugging in Eq. (Q) the appropriate luminosity functions 
(Granato et al. 2001). The total contribution of clustering 
to intensity fiuctuations has then been derived by adding up 
all the values of Cc{d) obtained for the different mass in- 
tervals and by also taking into account the cross-correlation 
terms between objects of different masses, according to the 
expression (where we have dropped the index C) 

3 

c^°^{0) = VcmcAe)- (8) 

The angular power spectrum of the intensity fluctuations 
can then be obtained as 

Ci^{W°\'^)^ / [5T{e)f PiicosO) sine de d(l), (9) 



with 



ST{e) = {(AT)')'''' = 



2 h 
hu - 2 



exp 



\ktTj 



( hv \ , f hv 



(10) 



which relates intensity and brightness temperature fluctua- 
tions. 

This analysis has been performed for the different cen- 
tral frequencies of the Planc k HFI, and Fi gure Is] shows the 
results in terms of 5Ti = ^Jl{} + l)Ci/27r (in units of K) 
at 2100 /xm, 1380 /xm, 550 /xm and 350 /xm. In each of the 
four panels of Fig. dotted curves are obtained by adopt- 
ing the 5(T detection limits of the different Planck channels 
{Slim = 300 mJy at 2100 /xm; Sum = 200 mJy at 1380 /xm; 
•S'lim = 450 mJy at 550 /xm; Sum ~ 700 mJy at 350 /xm) es- 
timated by Toffolatti et al. (1998). In Figure ^ we also plot 
predictions for temperature fluctuations at 850 /xm, where 
the contribution from star-forming spheroidals has been cor- 
rected for a mistake in MA2001. As in the former cases, dot- 
ted curves are obtained for a 5(t detection limit of the rel- 
evant Planck channel, corresponding to S'lim = 100 mJy, 
while dashed ones illustrate the case for ten times higher 
sensitivity {Sd = 0.1 Sum)- 

The solid curves in Figs. I and | show, for comparison, 
the power spectrum of primary Cosmic Microwave Back- 
ground (CMB) anisotropics corresponding to a flat ACDM 
cosmology, calculated with the CMBFAST code developed 
by Seljak & Zaldarriaga (1996). The relative importance of 
fluctuations due to clustering rapidly increases with decreas- 
ing wavelength. CMB anisotropies on small angular scales 
are exceeded at wavelengths A < 850/im. This high cluster- 
ing signal mostly comes from massive galaxies with bright 
fluxes, which lie at substantial redshifts and are therefore 
highly biased tracers of the underlying mass distribution. 
Also the negative K-corrections increase their contribution 
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to the effective volume emissivity [Eq. ([^] and therefore to 
the fluctuations. 

This implies that important information on the cluster- 
ing properties of faint sub-mm/far-IR galaxies (and hence on 
physical properties such as their mass and/or the amount of 
baryons involved in the star-formation process) will reside 
in the Planck maps at frequencies greater than 353 GHz 
where, however, the dominant signal is expected to come 
from interstellar dust emission. 

In order to quantify this last effect, we have calcu- 
lated the expected contribution of Galactic dust emission 
to background fluctuations. To this end, the dust emission 
maps constructed by Schlegel et al. (1998) by combining 
together IRAS and DIRBE data, have been extrapolated 
to all relevant frequencies assuming a grey-body spectrum 
with emissivity index m = 2 and average dust temperature 
Tdust = 18 K. 

Fig. 1^ shows the power spectrum of Galactic dust emis- 
sion averaged all over the sky (upper long-dashed curves). 
This signal dominates at all wavelengths due to contribution 
from the Galactic plane, highly contaminated by dust. Nev- 
ertheless, it is still possible to detect the clustering signal of 
sub- mm galaxies (and also recover the CMB power spectrum 
for A > 10^/im) if one restricts the analysis to high enough 
Galactic latitude regions (i.e. |fe| > 80°, lower long-dashed 
curves), which are least affected by dust emission. 

It turns out that, both for A = 850/im (see also Maglioc- 
chetti et al. 2001b) and A — 550/im, the dust contribution 
in these regions becomes less important than the one due 
to the clustering of unresolved sources at any I > 100. For 
A = 350/im the signal due to Galactic dust is instead found 
to hide all the other sources of background fluctuations at all 
Galactic latitudes, except possibly in the case of very high 
{I > 800) multipoles. 

Background fluctuations due to unresolved extragalac- 
tic sources have recently been measured by Lagache & Puget 
(2000) at 170 ^im in a held covered by the FIRBACK survey 
with ISOPHOT (Dole et al. 2001). 

As already discussed in this Section, the composition 
of counts at such relatively short wavelengths is substan- 
tially different from those considered so far because of the 
dominant presence of low-to-intermediate-redshift spiral and 
starburst galaxies. Calculations of the background fluctua- 
tions due to clustering of unresolved sources have then to 
take into account these two classes of sources, together with 
the population of spheroidal galaxies. 

We use again Eqs. (^) and ^ but, when dealing with 
spiral and starburst galaxies, an important difference in 
the form of the the bias factor focff (A/mm, 2) appearing in 
Eq. (^ has to be taken into account. In this case, in fact, 
feoff (Afmin, z) assumes the form: 

Kz) = l + ^ (11) 

(Fry 1996), independent of the mass of the haloes hosting 
such sources and merely a function of redshift - via the lin- 
ear growth rate D{z) - and of the local bias &o = b{z = 0) = 
o"8,9ai/o"8, where ag^gai is the rms fluctuation amplitude in a 
sphere of radius 8 h~^Mpc as measured in the local universe 
for each population. As a consequence, no distinction be- 
tween low-, intermediate- and high-mass objects is required. 
The choice for a bias evolution of the form given by 



Eq. ( [ll[ ) stems from a number of results (see e.g. Baugh 
et al. 1999; Maghocchetti et al. 1999; Magliocchetti et al. 
2000, to mention just a few) showing that bias models such 
as the one adopted for high-z forming spheroidal galaxies 
and illustrated by Eq. ^ do not provide a good flt for 2 < 
1, where most spirals and star-burst galaxies with 170 /im 
fluxes lying within the range relevant to this work, are found. 

We have estimated the amplitudes as,g from the clus- 
tering properties by spectral type determined by Loveday, 
Tresse & Maddox (1999) from Stromlo-APM survey data. 
We find as,g = 0.93 for galaxies with weak emission lines 
(spirals) and as.g — 0.66 for galaxies with strong emission 
lines (starbursts) . This leads to bo = 0.93 for spirals and to 
bo = 0.66 for starbursts, implying this latter population to 
be strongly anti-biased with respect to the underlying mass 
distribution. 

The contribution of clustering to C{9) has then been 
calculated for the three populations of spheroidals, spirals 
and starbursts. The angular power-spectrum A'^(fc), in the 
units used by Lagache & Puget (2000) (Jy^ sr^^), has been 
derived according to the expression (see Peacock 1997): 

poo 

A^(fc) = fcW c{e) Jo{ke) e de, (12) 
"'0 

where k is the angular wavenumber and Jo is the zero-th 
order Bessel function. The results for a detection limit Sd = 
135 mjy, corresponding to three times the confusion noise of 
the FIRBACK survey (Dole et al. 2001), are shown in Fig. 0, 
where the contribution of spheroidal galaxies is shown by the 
solid curves, while the dashed and dotted curves respectively 
show the predictions for spiral and starburst galaxies. 

As expected, the contribution of spheroidal galaxies at 
this wavelength is negligible when compared to those orig- 
inating from the clustering of both star-burst galaxies and 
spirals. Also, despite the strongly negative bias, the signal 
obtained for starbursts turns out to be stronger than the 
one predicted for spiral galaxies. This is due to the different 
evolutionary properties of the two populations, since no evo- 
lution is assumed for spiral galaxies, while strong evolution 
is instead allowed for in the case of starbursts (cf. beginning 
of this Section) . 

The short/long-dashed curve in Figure ^ represents 
the Poisson contribution Pp = 7400 Jy^/sr derived from 
the Guiderdoni et al. (1997) predictions, while the dashed- 
dotted curve illustrates the contribution to the observed 
power-spectrum of the Galactic cirrus (Lagache & Puget 
2000). According to our predictions, one then has that fluc- 
tuations stemming from the clustering of starburst galaxies 
are detectable above other sources of signal, and this hap- 
pens in the range of wavenumbers 0.04 < k (arcmin"^) < 
0.3. 

These findings, together with the predicted magnitude 
of the total (non-Poissonian) power spectrum for extragalac- 
tic sources (practically coincident with the dotted line in 
Figure |^ at 170 fim, are fully consistent with the measure- 
ments of Lagache & Puget (2000). We therefore conclude 
that an excess of signal (with respect to Poisson fiuctuations) 
due to clustering of unresolved/strongly evolving starburst 
galaxies, could possibly be discerned in the 170 jj,m FIR- 
BACK maps. 
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Figure 7. Predictions for the power-spectrum of intensity fluc- 
tuations due to clustering of galaxies fainter than a detection 
limit = 135 mjy at A = 170/im. The solid curves arc for 
the population of spheroidal galaxies (the higher one referring to 

-'^halo/A^sph = 100, the lower one to Mjj^io/Afjjpij = 10), while 
the dashed and dotted curves respectively show the predictions for 
spiral and star- burst galaxies. The short/long-dashed curve rep- 
resents a white noise power spectrum Pp = 7400 Jy^/sr, while 
the dashcd-dottcd one illustrates the cirrus confusion noise (both 
from Lagache & Puget, 2000). 



4 GRAVITATIONAL LENSING 

It is now well established (Franccschini et al. 1991, 1994, 
2001; Blain & Longair 1993; Pearson & Rowan- Robinson 
1996; Guiderdoni et al. 1997, 1998; Dwek et al. 1998; Blain 
et al. 1999b; Devriendt & Guiderdoni 2000; Pearson et al. 
2001; Takcuchi ct al. 2001; Rowan-Robinson 2001) that the 
coupling of the strongly negative K-correction at mm/sub- 
millimetre wavelengths - due to the steep increase with fre- 
quency of dust emission in galaxies - with the strong cos- 
mological evolution demonstrated by both ISO and SCUBA 
data (Elbaz ct al. 1999; Smail ot al. 1997; Hughes ct al. 
1998; Barger et al. 1999a; Bales et al. 2000) and by the in- 
tensity of the far-IR background (Puget et al. 1996; Fixsen 
et al. 1998; Schlegel et al. 1998; Hauser et al. 1998; Lagache 
et al. 1999, 2000), greatly emphasises high-redshift sources. 
This yields very steep counts which maximise the magni- 
fication bias (Peacock 1982; Turner et al. 1984) and lead 
to a high probability for such sources to be gravitationally 
lensed (Peacock 1982; Turner et al. 1984). As already dis- 
cussed by Blain (1996, 1997, 1998a,b, 1999, 2000), this cor- 
responds to a fraction of lensed sources expected to show up 
in the mm/sub-millimetre band which is much larger than 
what is found in surveys at other wavelengths. For instance, 
Blain (1998a) predicts about 0.6 to 5% of the point sources 
observed in the future by the High Frequency Instrument 
(HFI) on board of the ESA satellite Planck to be lensed. 

A distinctive feature of the astrophysically grounded 
model by Granato et al. (2001) is that both the lens- 
ing probability and the magnification bias at mm/sub- mm 
wavelengths are substantially higher than what other cur- 
rent, mostly phenomenological, models imply which also 
successfully account for SCUBA/MAMBO counts (Rowan- 
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Figure 8. Integral source counts per square degree for the 
model by Granato et al. (2001) at 550 (top), 850 (middle) and 
1380 /um. The light dotted lines, the heavy dashed lines, and the 
solid lines respectively refer to the unlensed, weakly lensed, and 
strongly lensed counts of forming spheroids only. The heavy dot- 
ted line shows the total counts (sum of weakly and strongly lensed 
sources), including contributions from populations whose counts 
are essentially unaffected by lensing (such as "flat" -spectrum ra- 
dio sources), as described in the text. 



Robinson 2001; Takcuchi et al. 2001; Pearson et al. 2001; 
Blain et al. 1999b). This is because, according to this model, 
most galaxies detected in blank-field SCUBA and MAMBO 
surveys are interpreted as massive ellipticals at a > 2, in 
the process of building their stellar populations with very 
high star-formation rates (typically from a few hundreds to 
~ IOOOMq yr~^). This leads to an extremely steep, essen- 
tially exponential, decline of the bright tail of mm/sub-mm 
counts of this galaxy population due to the fact that, in this 
model, dust emission from these objects rapidly fades away 
at z < 2 when the bulk of their star formation is essentially 
over (Cohen 2001). The bright counts therefore somehow 
reflect the high luminosity tail of the luminosity function 
which, in turn, reflects to some extent the fact that, in the 
hierarchical clustering scenario, massive halos are exponen- 
tially rare at high redshifts. Such a steep decrease of the 
number counts for fluxes lOmJy < Ssso^m ^ lOOmJy im- 
plies a large fraction of strongly lensed galaxies to appear at 
bright mm/sub-mm fluxes. 

We recall here the main aspects of our model for gravi- 
tational lensing, referring to our previous work (Perrotta ct 
al. 2001) for details. In the next Section we present quanti- 
tative predictions for the fraction of lensed spheroids to be 
compared with forthcoming observations. 

Lensing statistics is expressed by the probability for a 
source at redshift Zs to be lensed with magnification > /i: it 
is obtained by dividing the total lensing cross section by the 
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Figure 9. Integral source counts for the model by Rowan- 
Robinson (2001). The meaning of the lines is analogous to that in 
Fig. 8. In this case, however, the light dotted, the heavy dashed 
(superimposed to the dotted), and solid lines include the con- 
tributions of all the populations considered by Rowan-Robinson. 
The heavy dotted line also includes "fiat" -spectrum radio sources 
(see text). 



0.4 - 



0.3 - 



0.2 - 



0.1 - 



- 




Figure 10. Ratio of lensed to unlensed counts at 850 lira: the 
solid line refers to the model by Granato et al. (2001), the dotted 
line to the one by Rowan-Robinson (2001). 



area of the source sphere as 



dz 

dMa{p,z,z,,M)-^{z,M) 



(13) 



where r{z) is the comoving radial distance to redshift z, 
dV/dz is the proper volume element per unit redshift, and 
dn{z, M) /dM the comoving number density of the lenses. 
Since we are dealing with gravitational lensing by dark mat- 
ter haloes, we assume that the lens distribution follows the 
Sheth & Tormen (1999) mass function 



dn 
dM 



2aA^ po Scjz) 
TT M2 a{M) 



1 + 



\ y/aSc{z) 



2p 



diner 




dlnM 


exp ( 



(14) 



where the best-fit values of the parameters for our cosmo- 
logical model are a = 0.707, p = 0.3, and A ~ 0.3222. po 
is the mean mass density at a reference epoch to, which we 
assume to be the present time, and a is the variance of 
linear density fluctuations at the present epoch, smoothed 
with a spherical top-hat filter Wjj(fc) enclosing a mass M. 
5c (z) is the linear density contrast of an object virialising at 
z, linearly evolved to the present epoch. 

Quite independent of the lens model, P{p, z) decreases 
as for p ^ 1, hence the high magnification tail of 
the probability per unit magnification can be written as 
p{p, z) — —dP{n,z)/dp cx a{z) . Equation ( p^ as- 
sumes non-overlapping cross sections, which is satisfied in 
the P ^ 1 regime, i.e. when no more than a single clump 
causes lensing of a background source: this results in magni- 
fications markedly larger than 1, or strong lensing. Smaller 
magnifications, including de-magnifications, attain contribu- 
tions from the distribution of dark matter along the entire 
line of sight, resulting in weak lensing effects. The latter 
case can be represented by a Gaussian probability per unit 
magnification: 



p{p, z) = H{z) exp[-{p ~ p)^ /2al{z)] . 



(15) 



where the location of the peak, p, and the amplitude, H{z), 
are determined by the normalisation and flux conservation 
conditions obtained by integrating over all possible mag- 
nifications the combined (weak plus strong lensing) prob- 
ability distribution: J dpp{p,z) = J dppp{p,z) = 1. The 
transition between the weak and strong lensing regimes is 
set at a suitable magnification ^cut: a convenient choice is 
A*cut = 1 + 1.5(Tfi(2), yielding pcut ~ 1.5 — 2 for the redshift 
range of interest. 

The maximum magnification for strong lensing, pimax, 
determined by the intrinsic size of the source (Peacock 1982) 
and by the geometry of the lens configuration was found to 
lie in the range 10-30 in the case of our model (Perrotta et 
al. 2001). Here we adopt the conservative value pmax = 10. 

Lens density proflles have been modelled as Singular 
Isothermal Spheres (SIS) for masses M < 3.5 X 10" Mq, 
and with the Navarro, Frenk & White (1997; NEW) for- 
mula for larger masses. In fact, it is found that (Porciani 
& Madau 2000) such a "mixed" model provides a good fit 
to the observed statistics of the angular splitting of QSO 
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multiple gravitational images. However we wish to empha- 
sise that regarding statistical magnifications SIS and NFW 
profiles lead to comparable results (Perrotta et al. 2001), 
making the present conclusions very weakly dependent on 
the choice between SIS and NFW. 

Before going to the results, let us consider the magnifi- 
cation bias on a flux-limited source sample. The integrated 
source counts above a flux density threshold Si, of sources 
with a comoving luminosity function z) can be written 
as (e.g. De Zotti et al. 1996) 

N{S„)^ dz dL${L,z)r'^{z) -fsr-\ (16) 

where r is the comoving radial distance, and 

L^i„(^) = 47r(l + z) r\z) ^[(f ^^^^)^] ■ (17) 

The luminosity function modified by the magnification bias 
reads (e.g. Pei 1995): 

^'(L,-)-J^ d^^^<^>^^,z^ . (18) 

Lensing effects on the source counts are taken into account 
by replacing $'(1/, z) with <I>(L, z) in Eq. (p^). 



5 EFFECT OF GRAVITATIONAL LENSING 
ON SOURCE COUNTS 

We present here our estimates for the effect of lensing 
on integral number counts predicted using the model by 
Granato et al. (2001) and, for comparison, using the model 
by Rowan-Robinson (2001), taken as representative of the 
most successful phenomenological models. The model by 
Granato et al. (2001) includes three populations, namely 
spheroids, slowly-evolving starburst and evolving disk galax- 
ies (spirals); the model by Rowan- Robinson instead includes 
AGN, spirals, high-luminosity and low- luminosity starburst. 
In both models, we added the contribution due to "flat"- 
spectrum radio sources, whose counts (kindly provided by 
L. Toffolatti) are based on the model by Toffolatti et al. 
(1998). 

The contributions of different populations to the num- 
ber counts at fluxes probed by SCUBA predicted by the two 
models are markedly different: while, as discussed above, ac- 
cording to Granato et al. (2001) SCUBA sources are mostly 
massive dusty spheroids at high z, in the Rowan- Robinson 
(2001) model there is a large tail extending to \ow-z. Cor- 
respondingly, while the model by Granato et al. (2001) pre- 
dicts an essentially exponential decline of 850 fim counts 
with increasing flux density above a few mjy, the slope is 
substantially less steep in Rowan-Robinson's (2001) model. 

This different behaviour of the counts explains the 
vastly different fractions of lensed galaxies expected at 
bright fluxes, as quantified in Figs. ^ and ^ which show 
the integral number counts of lensed and unlensed galaxies, 
yielded by the two models, respectively, for the 550 (top), 
850 (middle) and 1380 /xm (bottom) Planck/HFI channels. 
The heavy solid lines correspond to strongly gravitationally 
lensed sources (magnifications fi > 2). Note that, for the 
model of Granato et al. (2001), forming spheroids give the 



relevant contribution to lensed counts, due to their expo- 
nential slope in number counts: thus, the solid lines in Fig. 
^ are given by the lensed spheroids, while starburst, spirals 
and radio sources are included only in the total counts (sum 
of weakly and strongly lensed sources), represented in Fig. 
^ by the heavy dotted lines. 

In the case of Rowan- Robinson (2001) all classes of 
sources (AGN, spirals and starburst) are taken into ac- 
count in the strongly- lensed counts of Fig. |^ (solid lines). 
The heavy dashed lines show the counts for these classes of 
sources, obtained by only allowing for weak lensing (magni- 
fications fj. < 2 and de- magnifications), while the thin dotted 
lines (almost coincident with the heavy dashed lines in the 
case of Rowan-Robinson 2001) show the counts computed if 
one ignores the effect of lensing. Finally, the heavy dotted 
lines give the sum of weakly and strongly lensed sources, 
also added with the contribution of radio sources, for which 
the effect of lensing is negligible. 

It is clear from Fig. |^ that, according to the model by 
Granato et al. (2001) and thanks to the strong magnifica- 
tions by gravitational lensing, even faint forming spheroids 
at substantial redshifts can be detected in relatively shal- 
low surveys, such as those provided by Planck/HFI. The 
quoted 5a detection limits for this instrument are 5*11111 = 
200,100,450 mJy at 1380,850,550 ^m, respectively. Note 
that, at 850 ^m, this implies that the expected strongly 
lensed objects at fluxes larger than ^um on the whole sky are 
~ 10^ in the Rowan Robinson scenario (2001) and ~ 4- 10"^ in 
the model by Granato et al. (2001). Strongly lensed form- 
ing spheroids may also be detectable by forthcoming bal- 
loon experiments operating at sub-mm wavelengths, such 
as BLAST (Devlin et al. 2000; www.hep.upenn.edu/blast/) 
and ELISA. Such surveys may then trace the large scale 
distribution of the peaks of the primordial density field. 

Note that the distribution of strongly lensed sources 
would extend to even higher fluxes in the case of magnifica- 
tions larger than our conservative limit ^max ~ 10 (magni- 
fications of up to 30 are possible according to the adopted 
model; cf. Perrotta et al. 2001). 

Fig. (jlo|) compares the fraction of lensed to unlensed 
sources as a function of 850 ^im fiux predicted by the 
two models considered above. According to the model by 
Granato et al. (2001), at 850 fim the ratio of strongly lensed 
to unlensed sources reaches the value of about 40% for fluxes 
slightly below 100 mJy, where the surface density of strongly 
lensed sources is of about 0.1 deg~^. On the other hand, the 
fraction of strongly lensed sources is always below 10% for 
the model of Rowan- Robinson (2001). It is then clear that 
gravitational lensing can help to discriminate different mod- 
els, all of them successfully reproducing the observed counts. 

In principle, clustered lenses can change the clustering 
properties of the sources in a flux-limited sample. The auto- 
correlation function of the lensing magnification pattern can 
be described as a convolution of the magnification pattern of 
a single lens with the correlation function of the lens centres. 
In our case, the magnification patterns of individual lenses 
have typical angular scales of at most a few arc seconds. The 
auto-correlation of lens centres, however, has a typical an- 
gular scale of arc minutes, and its amplitude is substantially 
lowered by projection because the lens population is taken 
from a large redshift range. Therefore, the effect of lensing 
on the clustering properties of sub-mm sources is negligible. 
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6 CONCLUSIONS 

A major challenge for current theories of galaxy forma- 
tion is to account for the evolutionary history of large 
spheroidal galaxies. Even the best semi-analytic models 
(Cole et al. 2000; Devriendt & Guiderdoni 2000) hinging 
upon the "standard" picture for structure formation in the 
framework of the hierarchical clustering paradigm, tuned to 
agree with detailed numerical simulations, are stubbornly 
unable to account for the sub-mm counts of galaxies as well 
for the (somewhat still controversial) optical evidence for 
most large ellipticals we see today were already in place at 
z ~ 2 and (almost) passively evolving since then (Daddi et 
al. 2000; Cohen 2001). These data are more consistent with 
the traditional approach whereby elliptical galaxies form 
most of their stars in a single gigantic star-burst followed 
by essentially passive evolution ("monolithic" models). This 
approach, however, is inadequate to the extent that it can- 
not be fitted in a consistent scenario for structure formation 
from primordial density fluctuations. 

A possible way out was proposed by Granato et al. 
(2001) who worked out a physically motivated scheme ac- 
cording to which the mm/sub-mm galaxies detected by 
SCUBA and MAMBO surveys are large proto-ellipticals in 
the process of forming most of their stars and of growing 
a quasar at their centres. In their model, the early evolu- 
tion of these galaxies and of quasars are tightly inter-related 
and feed-back effects (mostly due to supernova explosions) 
make large ellipticals to form as soon as the corresponding 
potential wells are in place, while the formation of smaller 
galaxies is delayed. Therefore, in this model, large ellipti- 
cals evolve essentially as predicted by "monolithic" models 
and the mm/sub-mm counts are successfully reproduced. 
The star-formation activity, powering the dust emission, 
quickly declines for z < 3 for the most luminous (most mas- 
sive) galaxies, while quasars reach their maximum luminos- 
ity. This naturally explains why very luminous quasars are 
more easily detected at mm/sub-mm wavelengths at red- 
shifts larger than that (~ 2.5) of maximum quasar activity 
(Omont et al. 2001). As indicated by the analysis of the 
latter authors, a substantial fraction of the observed far-IR 
emission is probably powered by the starburst in the host 
galaxy. 

In this paper predictions of the model are exten- 
sively investigated and compared with recent data. Atten- 
tion is focused, in particular, on two specific aspects. Since 
SCUBA/MAMBO galaxies are interpreted as massive galax- 
ies at 2 < 2 < 6, they are expected to be highly biased trac- 



ers of tint; maLLur disLribuLioii and Llierefory highly dusLyryd. 
The im plied angular correlation function is found to b e con- 
sistent with recent results by Scott et al. (2001), Peacock et 
al. (2000), Lagache & Puget (2000) as well as with the fluc- 
tuations in the SCUBA counts in different areas of the sky. 
Explicit estimates are presented for the power spectrum of 
temperature fluctuations due to clustering in Planck/HFI 
channels. 

A second specific prediction of the model is an essen- 
tially exponential decline of the counts at Sgsoum 10 mjy. 



Hints i 1 this direction that can possibly be discerned in 
recent SCUBA (Scott et al. 2UU1; Borys et al. 2UU1) and 
MAMBO (Carilli et al. 2000) surveys are noted. As an impli- 
cation of this prediction, together with the fact that sources 



are found at high redshifts, one has that both the gravita- 
tional lensing probability and the magnification bias on the 
counts are much higher than those derived for other current, 
phenomenological, models. We show that, according to this 
model, essentially all proto-spheroidal galaxies brighter than 
Sgsoum ^ 60-70 mJy are gravitationally lensed. Allowing for 
the other populations of sources contributing to the bright 
mm/sub-mm counts, we find that the fraction of gravitation- 
ally lensed sources may be ~ 40% at fluxes slightly below 
Sgso/jm = 100 mJy. If so, large area surveys such as those to 
be carried out by Planck/HFI or by forthcoming balloon 
experiment like BLAST and ELISA will probe the large scale 
distribution of the peaks of the primordial density field. 
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